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Abstract 

In the present paper, we give a series of definitions and properties of Lifting Dy¬ 
namical Systems (LDS) corresponding to the notion of deterministic diffusion. We 
present heuristic explanations of the mechanism of formation of deterministic diffusion 
in LDS and the anomalous deterministic diffusion in the case of transportation in long 
billiard channels with spatially periodic structures. The expressions for the coefficient 
of deterministic diffusion are obtained. 
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1 Introduction 

One-dimensional dynamical systems on the entire axis with discrete time are dehned by the 
recurrence relation 

Xn+l = f{Xn), (1) 

where / (a:) is a real function given on the entire axis and xq is a given initial value [H |T3j . 
Equation ([1]) determines a trajectory (xq, Xi,Xn, ■■■) = x in the dynamical system ([1]) 
according to the initial value Xq and the form of the function /. For the dynamical systems 
admitting the chaotic behavior of trajectories, the problem of construction of the entire 
trajectory or even of determination of the values of Xn for large n is quite complicated 
because, as a rule, the numerical calculations are performed with a certain accuracy and 
the dependence of the subsequent values of on the variations of the previous values is 
unstable. Moreover, from the physical point of view, the initial value Xq is specihed with 
a certain accuracy. Therefore for the investigation of the behavior of trajectories for large 
values of time, we can analyze not the evolution of system ([T]) but the evolution of measures 
on the axis generated by this evolution. 

If a probability measure /tq (a normalized measure for which the measure of the entire 
axis is equal to 1) with density po ^ ho(^) = / Po{.x) dx is given at the initial time, then, for 

A 

a unit of time, system ([T]) maps this measure into pi : pflA) = p,(/“^(A)), where f~^{A) is 
a complete preimage of the set A under the map /. The operator mapping the measure po 
into the measure pi is called a Perron-Frobenius operator. The Perron-Frobenius operator 
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T is linear even for nonlinear dynamical systems ([T]) and maps the density po of the initial 
measure into the density p\ of the measure pi as the integral operator with singular kernel 
containing a Dirac 5 function: 


pl{x) = J^po{x) 


- f{y))pQ{y) dy 


E 


Vk&f 


\f'{yk)\ 


PoiVk 


( 2 ) 


The investigation of the asymptotic behavior of the density p„ = as n —)■ oo is 

reduced to the investigation of the behavior of the semigroup There are examples of dy¬ 
namical systems CQi with locally stretching maps / for which the densities pn are asymptoti¬ 
cally Gaussian as n —)• oo independently of the choice of the density of the initial probability 
measure. In this case, it is said that deterministic diffusion occurs in the dynamical system 
©• 

The aim of the present paper is to consider examples of dynamical systems with deter¬ 
ministic diffusion. We restrict ourselves to the so-called lifting dynamical systems (LDS) 
with piecewise linear functions /(x) in DS ([T]). 

We briefly consider the mechanisms of appearance of the anomalous deterministic diffu¬ 
sion in the process of transportation in long billiard channels with spatially periodic struc¬ 
tures. 


2 Lifting Dynamical System 


Consider the dynamical system ([T]) on the entire axis, where the function f{x) is given on 
the main interval Jq = [—|, |). The function s(x) = f{x) — x with x G Jq has the sense of a 
shift of the point x under the map /. We split the entire axis (— cxd, oo) into disjoint intervals 
h = [k — + ^), where k E Z are integers. Assume that the function / is extended from 

the interval Jq onto all intervals Ik so that the shift of a point under the map / on each 
interval k is the same as on Jq, that is s{k -|- x) = f(k + x) — {k + x) = s(x) = /(x) — x. This 
gives a periodic shift function s(x) and implies the periodicity of the function / with lift 1 

f{k + x) = k + f{x), \x\ k e Z. (3) 


The dynamical system ([T]) with a function / satisfying property (jS]) is called a LDS. This 
dynamical system (DS) is well known and thoroughly investigated [U El IS]. 


Lemma. Let a function f defined on the interval Iq = [— 2 , 2 ) ® piecewise monotone 

m 

stretching function on a finite partition of the interval /o = U ^oj, i-e., there exists 


A > 1 such that 


j=i 


\f{x) - f{y)\ > A|x - y\ 


( 4 ) 


for x,y G lo^k, (ind, in addition, the function /(x) has finite nonzero discontinuities at the 
points of joint of the intervals Iqj. Then the periodic extension ^ with lift 1 specifies a 
locally stretching function f on the entire axis. 


Proof. If the points x and y belong to the same interval /o,fc, then, by virtue of dlj). 


\f{x) - f{y)\ = |/(x - k) - f{y - k)\> A|x - y\. 
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If the points x and y belong to neighboring intervals, then the fnnction has a discontinnity 
at the point of joint of these intervals and the valne of this discontinnity is not smaller than 
a certain a > 0. In view of the monotonicity of the fnnction /, there exists £o > 0 snch that, 
for \x — y\ < the points x and y belong either to the same interval /o,fc or to neighboring 
intervals. Hence, 

\f{x)-f{y)\>l>Xe, 

where £ = min(£o, ^)- Thns, for \x — y\ < e, we have \ f{x) — f{y)\ > A|a; — y\. □ 

Let X = {xq, Xi,..., Xn, •••) be a trajectory for the LDS ([I])-(l3]). By fc, we denote the nnmber 
of the interval Ik = [fc — containing a nnmber a G /a, as follows: k = [a) is the nearest 

integer for the nnmber a. Then the trajectory x is associated with an integer-valned seqnence: 
M = ([a:o), [ 3 ^ 1 ), •••, [xn)) ■■■) = {^ 1 , 1112 , ■■■) of the nnmbers of intervals containing the 

valnes Xn- The seqnence M is called the ronte of the trajectory x and contains the nnmbers 
of intervals snccessively “visited” by the phase point in the LDS [8]. 

Proposition 1. Let the fnnction / determining the LDS ([I])-© be stretching. Then the 
trajectory of the LDS is nniqnely determined by its ronte. 

Proof. Assnme that two initial conditions and a;® generate trajectories of the LDS with 
the same ronte. This means that the points Xn'^ and Xn^ he in the same interval Im„- In view 
of the stretching property of the map /, the ineqnality |a:^+i — > A|a;n^ — a;n^|, A > 1, 

is trne. Since any two points from the same interval differ by at most 1, the snccessive 
application of this ineqnality yields the ineqnality |a;o^^ — a;g^^| < ^ for any n. Since A > 1, 
this implies that a;o^^ = x^q\ □ 

Definition 1. If any integer-valued sequence in the LDS is a route of a certain trajectory 
and the trajectory is uniquely determined by its route, then we say that the LDS possesses 
the Bernoulli property [8]. 

Proposition 2. If a function / defined in the interval Jg = [—|, |) is odd, continuous, 
monotonically increasing, stretching, and unbounded, then the LDS ([I])-® corresponding 
to this function / possesses the Bernoulli property. 

Proof. The uniqueness of reconstruction of a trajectory according to its route is proved in 
Appendix [U Let M = (mi, m 2 ,..., rUn,...) be an arbitrary integer-valued sequence. In the 
interval Imi , we consider a sequence of embedded contracting closed intervals f/„ constructed 
as follows: 

Let be the closure of the preimage of the interval under the map / lying on Imn-i 
: = f~^{Im„)- Since the map / is stretching, the length mesf/*-^^ < A”^. Let be the 

preimage of in the interval /m „_2 and let, by induction, be the preimage of in 

the interval Im„_k- The sets Un = C Imi and raesUn < A”. As n —)■ 00, we get a system 
of embedded closed sets Un+i C Un in the interval Imk and mesUn —)■ 0 as n ^ cx). The 
common limit point Xq for all Un lies in and its trajectory x = (xg, f{xo ),..., /""(xg),...) 
has the route M. □ 
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3 Markov Partition of the Phase Space of LDS. Piece- 
wise Linear LDS 

The system of intervals Ik = [k — k + ^), k E Z, for the LDS ([I])-© can be regarded as 
a Markov partition of the phase space [5]. Consider a finer Markov subpartition. Assume 
that the main interval Jq is split into finitely many m subintervals Iqj = where 

xq = —\-iXj < Xj+ijXm = = 0,1,2, and Joj = [xj-i^Xj). An integer-valued 

shift of this partition leads to the decomposition of the intervals Ik into subintervals Ik,j = 
[k -|- Xj-i, k + Xj). Hence, the entire axis is split into intervals {hj}, k E Z, j = I, 2, m. 

Definition 2. We say that the Markov partition {Ik,j}k&z,i<j<m of the entire axis is consis¬ 
tent with the LDS ([I])-(l3j) if the LDS maps, for a unit of time, any probability measure with 
constant densities in each set Ik,j into a measure with constant densities in each interval Ikj- 

Definition 3. We say that a function / specifying the LDS ([I])-® is consistent with the 
Markov partition {Ik,j}kez,i<j<m dehned with the help of the numbers —2 = xq, Xi, ...,Xm = 
2 if it is linear and nonconstant on each interval Ikj and, at each end of the interval Ikj, 
takes values equal to an integer plus one of the numbers Xj, j = 0,1,..., m. 

Example 1. Let f{x) = Ax be a linear function in the interval Iq = [—|, |), where A = 2/-I-1 
is an odd number. Then this function is consistent with the Markov partition {Ik}k£Z, Ik = 
[k — k + I). 

Example 2. Let f{x) = Ax, where A = 2/ is an even number. Then this function is 
consistent with the Markov partition {Ik,±}kez, where Ik,+ = [k,k + \) and Ik - = [k — \, k). 

An important example of the function / satisfying Dehnition [3] is given by the following 
assertion: 

Proposition 3. Assume that a function f{x) determining the LDS in the interval Jq is 
piecewise linear, takes half-integer values at the ends of each its linear pieces, and moreover, 
f'{x) 7 ^ 0 almost everywhere. Then, by Dehnition [31 the function / is consistent with the 
Markov partition {Ik}\k&z- 

Proof. The proof follows from the verihcation of the conditions of Dehnition [3] for points 
{xj}, that is for all points at which the function / takes half-integer values. □ 

A piecewise linear function / from Proposition [3] will be called a piecewise linear function 
taking half-integer values at the ends of the linear pieces. 

The equivalence of Dehnitions |2] and [3] for LDS yields the following statement: 

Proposition 4. For the Markov partition {In,j}nGZ,\j\<m to be consistent with the action of 
the LDS ([I])-(|3]) by Dehnition [21 it is necessary and sufficient that the function / determining 
the LDS ([I])-(13]) be consistent with this Markov partition by Dehnition [31 

Proof. We now show that if / satishes the conditions of Dehnition [31 then the LDS transforms 
the probability measure with constant densities on j into a measure with constant densities 
on In^ki he., the LDS satishes the conditions of Dehnition [21 Since operator ((21) is linear, it 
suffices to consider the case where the density of the initial measure is constant on the hxed 
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interval InoM- interval, the fnnction f{x) is linear and nonconstant and maps InoM 

into the nnion IJ ^ of several neighboring intervals. Therefore, the inverse map of these 
intervals is linear and, hence, the measnre has a constant density in each of these intervals. 

Let the conditions of Dehnition [2] be satished. If the density of the initial measnre is 
constant in the interval /no.fco) then the map / gives a new measure /ii on the entire axis. It 
is clear that supp/r coincides with the closure of a certain union f~^{In,k) = U-^oy because, 
otherwise, there exist an interval and its part A such that ^ 0 and 

fii{A) = 0. This contradicts the condition of Dehnition [2] according to which the measure 
fii has a constant density. Let 1^ ^ be one of the intervals /(/no,fco) let B = 

Since B C Ino,ko, the initial measure on B has a constant density, / bijectively maps B 
onto Injti and, according to the Perron-Frobenius operator, the density pi{x) of the measure 
Pi on is expressed via the density po of the initial measure on Ino,ko by the equality 
Pi{x) = . Thus, we arrive at the conclusion that f'{x) = const on the interval B and 

that the function maps the ends of the interval B into the ends of the interval I^k- Thus, 
the function / satishes the conditions of Dehnition [31 □ 

The following question arises: Is it possible to construct Markov partitions of the en¬ 
tire axis consistent with the linear function f{x) = Ax for the values of A other than in 
Examples [H and [2]? 

Since the linear function f{x) = Ax is odd, the numbers Xj specifying a consistent Markov 
subpartition of the interval Jq are symmetric about the middle of the interval Jq. Hence, it 
is sufficient to dehne solely the positive values of Xj and indicate that the number m of 
subintervals {Io,j}i<j<m is even or odd, because s = 0 belongs to the set {s^} for even m 
and does not belong to this set for odd m. 

We enumerate the ends of the intervals {Ikj} located on the positive half axis in the 
order of increase starting from the hrst positive number. This yields a sequence 0 < sq < 
Si < S 2 < ... < Sk < .... Since the Markov subpartition {hj} is consistent with the LDS 
with the linear function f{x) = Ax, we get 

ASi Si_|_^j, As 2 Asi_|_^j_|_^2 ) ■ • •) As^ '^l-|-ni-|-n2-|-...-|-n,^ • (b) 

Here, rh is expressed via m of the form m = 2rh for even m and of the form m = 2m—1 for odd 
m. Hence, each Markov subpartition of the axis consistent with the linear function f{x) = Ax 
specifying the LDS ([I])-([3]) is associated with two parameters: the parity of the number m of 
subintervals {Iqj} (i.e., s = (—1)™) and the integer-valued vector n = (ni,n 2 , ...,nm). These 
quantities (s,n) are called the parameters of the Markov subpartition {hj}- 

Proposition 5. The parameters of the Markov subpartition {Ik,j}, he., the parity of the 
number m and the integer-valued vector n = {ni,n 2 , uniquely dehne the quantity A, 

the slope of the linear function f{x) = Ax, consistent with the Markov subpartition {/fej} 
and a collection of numbers specifying the ends of the intervals Ikj = [k + Sj, k + Sj^i) 

of this subpartition. 

Proof. Let the number of components of the vector n be equal to m. Then the number m 
of subintervals {/o,j} in the Markov subpartition is given by the equality m = 2m for even 
m and the equality m = 2m — 1 for odd m (the parity of m is specihed by the quantity 
s = (—1)™). We always have Sm = |. For odd m, the sequence si, S 2 , ■■■, Sk,... has the 
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form Si, S 2 , Sm+l, 1 ~ Sm-l, 1 ~ Sm- 2 , 1 — Si, 1 + Si, 1 + S 2 , ■ ■ •, 1 + S^-l, |, 2 — S^-l, 

and is explicitly expressed via Si,S 2 , ■••,Sm = Hence, we can explicitly express each term 
of the seqnence Si, S 2 ,..., s^,... in terms of an integer and one of the numbers Si, S 2 ,..., s^. 
Thus, equalities ([5]) turn into a linear system for si,S 2 ,...,Sm whose coefficients are either 
integers or the quantity A. The consistency condition for this system can be formulated as 
the equality of the determinant of this system to zero. This gives the following algebraic 
equation for A: 

Rn{A) = 0 . 

If we hnd A, then all si,S 2 ,...,Sm can be uniquely determined from ([5]). Similarly, we 
consider the case of even m. In this case, the sequence si, S 2 ,..., Sk ,... contains integer values 
and has the form si, S2,..., Sm = ^, 1 - s^-i, 1 - Sm-2, ..., 1 - si, 0,1 + si, 1 + S2, .... □ 

An implementation of this scheme by a constructive example is given in the next Propo¬ 
sition. 


Proposition 6. A linear function f{x) = Ax determining the LDS ([IJ-Q will be consistent 
with the Markov partition {[fc —^ + 0) + ^ + |)}A:ez, if for integers 0 < m < n 
and the values ei = ±1, £2 = ±1 we have 


2 n + 82 + \/ (2n - £ 2 )^ + Sm^i 
“ 2 ’ 

^ 2m 

2n - £2 + \/ (2n - £ 2 )^ + Smei 

Proof. This is a consequence of the equalities ([5]), in this case having the form 


Af,=m + £ 2 ^, 

A 

- = n + eif. 

□ 


It is worth noting that the set of values of the slope A of the linear function f{x) = 
Ax,x G Jo = [—|, |), for which it is possible to construct Markov partitions consistent with 
the LDS ([I])-® is everywhere dense on the half axis (2, cxo) 12]. 


4 Deterministic Diffusion 

We consider the LDS ([I])-© for which the collection of intervals {Ik}kez, h = [k — k + ^), 
forms a Markov partition of the phase space consistent with the action of the LDS. This 
means that, for a time unit, the LDS maps the measure po with unit density in the interval Jq 
into a probability measure with constant densities > 0 in the intervals R and ^ = 1. 

k&Z 

As a result of multiple application of the LDS, the initial measure po is transformed into a 
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measure with constant densities in the intervals Ik- Let Pk{n) be the density of the measure 
in the interval Ik after the n-fold action of the LDS upon the initial measure fiQ. Then 

Pk{n + 1) = ^ pk-iPi{n) (6) 

i&z 

and Pfc(O) = 5kfl due to the choice of the initial measure /iq. The asymptotic behavior of 
the quantities Pkiji) for large n, as solutions of Eq. ([6]), is described by the well-known 
central limit theorem of the probability theory [SI E]. Indeed, if we consider the sum ^ = 
+ ^2 + ••• P oi n independent random variables each of which takes only integer values 
k with probability then we get Eq. ([6]), where Pk{n) is the probability of the event that 
^ takes the value k. 

Theorem 1. Let the numbers {pk}k&z take nonnegative values such that pk = 1, let the 

k&Z 

greatest common divisor of the numbers k with pk > 0 be egual to 1, and let there exist the 
first and second moments 


(Ti = ^ kpk, = X] (7) 

k£Z k£Z 

Then the solution of Eg. & as n —)■ oo has an asymptotics 

Pfc(n)-2-"" -)■ 0 (8) 

ay ZTin 


uniformly in k. The character of convergence in m depends on additional conditions, e.g., 
on the existence of the third moment (the Laplace theorem). 


Proof. We now briefly present the well-known scheme of the proof of Theorem [T] based on 
the fact that Eq. (E]) is a difference analog of the convolution equation. As a result of the 
Fourier transformation, the convolution turns into the product. Let P{X,n) = ^ Pk{n)e^^^ 

k&Z 

be the characteristic function of the solution of Eq. ([H]), P(A) = ^ pkC^^^. Thus, we get the 

k&Z 


following formula from Eq. (E]): 


F(A,n) = [P(A)]"P(A,0). 


(9) 


Since the initial measure is concentrated on Jq and its density is constant, we conclude 
that P(A, 0) = 1 and 

TT 

PfcH = ^ j[P{X)Te-^^^d\ (10) 

— TV 

As A ^ 0, P(A) = X] PkP^^ = 1-1- iaiX — + o(l) = -|- o(l). Hence, 

k&Z 

[P{X)Y = + o{l). 

Substituting this result in f[T0|) and integrating along the entire axis, we obtain ([8]). 

□ 
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The degree of closeness of two probability measures /i and v with densities 'p{x) and q{x) 
on the axis is often estimated by the value d{fi, v) of the deviation (uniform in x) of the 
distribution functions 

X 

P(x) = j p{s)ds = fx{{-oo,x)), 

— OO 
X 

Q{x) = J q{s) ds = u{{-oo,x)), 

— OO 

i.e., 

d{fi, v) = sup l-P(x) — Q{x)\ (11) 

X 

Definition 4. We say that two sequences of measures pn and z/„ are asymptotically equiv¬ 
alent as n —)■ cxD if lim d{nn, Vn) = 0. For normal measures Vn with variances and means 

n^oo 

^n, i-G-; in the case where the density of measures has the form of a Gauss curve 

1 

Qnix) = — , ( 12 ) 

we say that the sequence of measures pn is asymptotically normal as n —)■ oo. In addition, if 
= Dn, then we say that the sequence of measures /i„ determines a normal diffusion with 
diffusion coefficient D. If the variance depends nonlinearly on n, then the deterministic 
diffusion is anomalous. 

The result of Theorem [T] for the LDS can be interpreted as follows: The initial measure 
with unit density in the interval Jq can be regarded as randomly specihed initial data Xq for 
the LDS uniformly distributed over the interval Jq. Then, for large times, as n ^ oo, the 
position of Xn is randomly distributed according to the normal law with mean value and 
variance cr^. 

In other words, in this case, we have the deterministic diffusion with the diffusion coeffi- 
2 

cient D = ^. Our hrst aim is to study this phenomenon. 

First, we reformulate the result of Theorem [1] for the LDS ([I])-([3]) with a piecewise-linear 
function /. 

Theorem 2. Assume that a function f determining the LDS ([I])-([3) is a piecewise linear 
function taking different half-integer values at the ends of all linear parts. 

Suppose that there exists 


1 

2 

D=\ J \nO\''dx-L (13) 

_ 1 
2 

Then, after n iterations in the LDS, the initial measure po with unit density in the 
interval Iq is asymptotically mapped into a measure with normal distribution and the diffusion 
coefficient D. 





Proof. In view of Theorem [T] and Propositions |3] and IH it is necessary to show that 

1 

2 

= / \f{x)\‘^dx- (14) 

k&Z \ 

2 

This can readily be proved becanse the integral of the piecewise linear fnnction / can 
easily be taken. If the fnnction /(x) is linear in the segment [a, h] and takes values | 

b 

and fc + I at the ends of this segment, then J \ f{x)\^ dx = {hf + — a)- 

a 

We get pk = pi{Ik) = p{f~^{dk))- The set f~^{Ik) fl-^o is formed by several subintervals 
Iqj of the interval Jq in which the function /(x) takes values from the interval [j — j + ^)- 
Hence, / \f lx)\^ dx = {k'^ + polloj) = + ^)Pk- 

f-Hik)r\h 

As a result of summation over fc, we obtain 


1 

2 


1 

2 


\f{x)\^dx = '^{P + ^)pk, 
kez 


which is equivalent to flTT)) . □ 

Example 3. Let /(x) = Ax, where A is a positive number. Then function / satishes the 
condition of Theorem [2] if and only if A is an odd number and A > 1. Formula fll3p gives 
D = 1/24(A^ — 1), which of course, coincides with the formula known in the literature (see, 
for example, (23.21) in |3] and references therein.) 

Example 4. (zig-zag map) On the interval [—1/2, 1/2], let us consider an odd piecewise 
linear function / that takes the half-integer value /(^) = p-|-l/2ata point ^ (0 < ^ < 1/2). 
Let /(O) = 0 and /(1/2) = 1/2. Then the diffusion coefficient, according to f[T3ll . has the 
following form: 

D = ^{2p+l-20. 

It argees with known results (see na and the references therein.) 

By using the results of Theorems [1] and [2], we can give the following definition of deter¬ 
ministic diffusion for the dynamical system (HJ: 

Definition 5. We say that the one-dimensional DS ([T]) has a deterministic diffusion if, for 
any initial probability measure po with bounded density, there exist a sequence of numbers 
cr^ > 0 and such that the sequence of measures /r„ = F"’Pq obtained from the initial 
measure by the n-fold action of the DS is asymptotically equivalent, as n —)■ oo, to a sequence 
of normal measures with variances and mean values 

The main problems connected with the deterministic diffusion for DS is to establish the 
fact of existence of this diffusion in terms of the functions / specifying the DS ([T]). The 
construction of an efficient algorithm for the determination the coefficient of deterministic 
diffusion D and drift seems to be an important problem, especially from the viewpoint of 
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applications. At present, a series of expressions is dednced and various numerical methods 
for the analysis of the dependence of the diffusion coefficient on the form of the functions / 
are developed. Note that the dependence of the coefficient of deterministic diffusion for the 
LDS ([I])-® on A is quite complicated (nowhere differentiable fractal dependence) even for 
the linear function /(x) = Ax in the interval Iq = [—^, |), [HIS]- 

We now present heuristic arguments for the existence of deterministic diffusion for the 
LDS with the linear function /(x) = Ax, A > 2, in the interval Iq = [—|, |). In this case, 
the LDS ([I])-([2]) can be represented in the form 


x„+i = x„ + (A - l){x„), (15) 

where {x„) = x — [x) is the fractional part of the number x and [x) is the nearest integer for 
the number x. Equation fll5p can be represented in the equivalent form 

n 

x„+i = Xo + ^(A - l){xfc). (16) 

k=l 

If the map / is stretching, i.e., A > 2, and the initial value Xq takes values in the interval 

Iq with a certain (e.g., constant) probability density, then we can assume that the quantities 

{xfc), the fractional parts of Xk are uniformly distributed over the interval Iq and independent 

for different k. The rigorous substantiation of the uniformity of distributions of the fractional 

parts for different stretching maps can be found in [7]. According to flTB]) . the quantities x„ 

can be regarded as the sum of Xq and n independent identically distributed random variables. 

Thus, by the central limit theorem, the quantities x^, as n ^ cx), are distributed according 

to the normal law with zero mean value and the variance equal to the sum of variances of 

the terms. Since the variance of {x), regarded as a variable uniformly distributed in the 

1 

interval [—|, |), is equal to cr^ = / x^ dx = we get ct^(x„,+i) = . This yields the 

_ 1 
2 

approximate relation for the diffusions coefficients D = for any linear map /(x) = Ax 

in the LDS ([I])-([3]). 

Let us turn to the case where the Markov subpartition is consistent with the LDS ac¬ 
cording to De£nition[21 In this case, by Pkj{n), we denote the densities in each subinterval 
Ik,j of the interval Ik, k ^ Z, j = 1,2, for n iterations. If we consider the collection 

of Pk,j{n) j = 1,2, ...,m, as the components of the vector Pk{n) = col {Pj^^ then 
we get an analog of Eq. ([6]), where the quantities Pk{n) are vectors, pk = {Pk,i,j}\i^j=i is the 
transition matrix for the Perron-Frobenius operator, and Pk,i,j is the density of measure in 
Ik^i in the case of single application of the LDS to the measure with unit density on Ig j. 
This vector analog of Eq. ([6]) is also well studied and the solution Pk{n) leads to the normal 
distribution in the m-dimensional space 013]. 

Theorem 3. Assume that the Markov subpartition {Ik,j}kez,i<j<n of the axis is consistent 
with the action of the LDS ([I])-([3) o^nd that the function /(x) is stretching and maps the 
interval Iq = [—1,|) into a finite interval [a,b) of length greater than 2. Then, after the 
n-fold action of the LDS, the initial measure Pq in Iq with bounded density is transformed 
into a measure pn asymptotically eguivalent, as n ^ oo, to a normal measure with densities 
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Pk,j{,n) in the intervals Ikj 


Pkji'n) 


aj 

-g 4Dn 

2\/ ttDu 


j = 1,2, 


(17) 


where is the drift, D is the coefficient of deterministic diffusion, and the parameters 
aj >0, j = l,2,...,m, specify the distributions of densities in the subintervals Ikj of the 
intervals Ik,k E Z. 


Proof. As already indicated, the vectors Pk{ri) = col {Pk{ri, 1),..., Pk,mijk)) satisfy the equa¬ 
tion 


Pkin + 1) = Y^ PjPk-jin), (18) 

where the matrix Pj is expressed via the translation matrices in the considered LDS. Note 
that the matrix E = 'Yhj Pj is equivalent to the stochastic irreducible matrix DED~^, where 
is a -D-diagonal matrix with the lengths of subintervals Joq, on the diagonal. As a 

result of the Fourier transformation, Eq. flT8|l is transformed into the difference equation 

P{X,n + l) = P{X)P{X,n), (19) 

where the matrix P{X) = ^ PjC^P and the vector P{X,n) = ^ Pk{n)e^^^. 

j k 

The solution of Eq. flT^ is explicitly expressed via the eigenvalues and eigenvectors 
(including the adjoined vectors in the case of multiple eigenvalues) of the matrix P{X). If 
2 ;(A) is the maximum eigenvalue of the matrix P{X) (this eigenvalue is simple), then the 
solution of Eq. flT^ as n —)■ oo can be represented in the form 

F(A, n) = 2 :"(A)a(A) -|- o(l), n ^ oo, (20) 


where a(A) is the eigenvector of the matrix P{X) corresponding to the eigenvalue ^(A). Note 
that ;2(0) = 1 and all components of the vector 0(0) are positive. If we perform the inverse 
Fourier transformation, then we get relation flTT)) with D = —|^ 2 ;(A)|a =05 ^ = i^^^^ 
from fl20|) . 


ax 


A =0 

□ 


Remark 1. By using Theorem [3l one can deduce the explicit “parametric” dependence of 
F on A for the linear function /(x) = Ax specifying the LDS. In this case, the role of 
parameters is played by the characteristics (s, n) of the Markov subpartition {Ik,j} consistent 
with the LDS, that is by the parity of m and the integer-valued vector n. By using these 
parameters, we can explicitly construct three polynomials Rn{x), Pn{x), and Qn(x) with 
integer-valued coefficients such that the slope A is the maximum root of the polynomial 

Rn{x), i.e., Rn{x) = 0 (see Proposition [5]) and D = 


P^ 


X 


Qf 


[X 


Example 5. 


As an example, we consider the case of LDS with a linear function /(x) = Ax, where 
A = 2s is an even number. In this case, the Markov subpartition Ik is formed by two 
subintervals Ik,+ = [k,k + |) and R- = [k — |,fc). If the hrst components of the vectors 
are referred to the intervals 1^,+ and the second components are referred to the intervals 
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0 1 
0 1 


Ik-, then Po = 


A 


1 0 
0 1 


Pj 


1 

A 


Ps = 


0 0 
1 0 


p-j = 


A 


, and 


A 


0 1 
0 0 


, j = 1,2, ...,s — 1. In this case, the characteristic function P{X,n) = 


EPk jg vector-valued and 


pw = E p., 


sinfs 




0 2 


2s sin I 


e 2 


The determinant of the matrix P{X) is equal to zero and the trace spp(A) =-^ cos —(s — 

s sin j 2 

1). Hence, the nontrivial eigenvalue z(A) of the matrix P{X) coincides with the trace of the 
matrix P{X). This yields the following well-known explicit expression |3] for the diffusion 
coefficients for even A = 2s: 


D 


1 

2 ^ 


spp(A)| 


A=0 — 


(A-l)(A-2) 

24 


( 21 ) 


As it follows from the relations of Example [3] and fl 2 T]) for the coefficient of deterministic 
diffusion in the case of a linear function f{x) = Ax with integer A, -D(A) is a monotonic 
function of A. As A increases, i.e., the degree of stretching of the map f{x) increases, the 
coefficient of deterministic diffusion D{A) increases. However, if we compare the relations 
for even and odd A, then, e.g., we get Us = A for A = 3 and 114 = A for A = 4. Thus, at 
hrst sight, these results seem to be intuitively strange. 

We consider this problem in more detail. Assume that the initial probability measure 
is uniformly distributed over the interval [|,2), that is its density is constant and equal to 
2 . As a result of the one-time action of the LDS, we get a measure with constant density 
in the interval [A, 2) for the mapping with A = 3 and, in the interval [0, 2) for the mapping 
with A = 4. Clearly, the variance of a measure with constant density in the interval [0, 2) 
is greater than the variance of a measure with constant density in the interval [^, 2 ), which 
agrees with our assumption that the higher the degree of stretching, the greater the variance. 
A similar picture is observed in the case where the initial measure has a constant density in 
the interval [—2, —|). For A = 3, the mapping realized by the LDS leads to a measure with 
constant density in the interval [— 2 , — A) and the mapping with A = 4 yields a measure with 
constant density in the interval [—2,0). 

However, if the initial probability measure has a constant density in the union of intervals 
[— 2 , —|) 1J[|, 2 ), then, in view of the linearity of the Perron-Frobenius operator, for A = 3, 
we get a measure in the union of intervals = [— 2 , —A) IJ[A^ 2 ], whereas for the mapping 
with A = 4, the measure has a constant density in the interval A 4 = [—2,2]. Clearly, the 
variance of the probability measure on A^ is greater than the variance on ^ 4 . Thus, the 
higher degree of stretching not always corresponds to a greater variance. If stretching is 
directed toward the mean valne, then the variance decreases. 

The relations for the values of D for even and odd A can be replaced by a single relation 
valid for all A by introducing a function a;(A), namely. 


D(A) = i(A-l)(A-i.,(A)), 


( 22 ) 
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where the function a;(A) is equal to 2 for even A and to -1 for odd A. For any A, the function 
a;(A) is characterized by a complex fractal dependence on A caused by the fractal dependence 
of D on A (see Hi). The function a; (A) can be approximately regarded as a 2-periodic 
function such that Ci;(A) = 2 — 3|A — 4| for A G [3, 5]. The plot of this function is depicted in 
Fig. 1 and gives the first approximation of the behavior of D as a function of A. 



Figure 1: Graph of the function D{A). 

Let us consider some examples of a linear function f{x) = Ax, that illustrated the results 
of theorem [31 

Example 6. 1. Markov subpartition of m = 3 subintervals: = [k — ^,k — ^),Ik ,2 = 

[k — k + ^), = [k + ^, k + ^), k ^ Z. 

2. A system of two equations for ^ and the value of slope A : 

3. An equation for A : A^ — 4A + 1 = 0. 

4. A value of A : A = 2 + \/3 ~ 3.73. 

5. A value of ^ ^ ^ 0.13. 

6 . A matrix P(A) : 

/ 1 e*"'' + \ 

F(A) = - 1 

y g-*^ e-2*A ^ giA J 

7. A value of greatest eigenvalue 2 ;(A) of matrix P(A) : 

2 ;(A) = + cos A + \/cos^ A + 2 cos A). 
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8. A value of deterministic diffusion coefficient: D = ^ ^ 0.29. 

D 

9. Values of densities aj on subintervals Ikj : ai = as = 7 ^ 1.08, 0:2 = ~ 0.79. 

Example 7. 1. Markov subpartition of m = 3 subintervals: 

7fc,i = [k-\,k- 0,4,2 = [k - i,k + = [k + i,k + \),k e Z. 

2. A system of two equations for ^ and the value of slope A : 

3. An equation for A : A^ — 6A + 3 = 0. 

4. A value of A : A = 3 + \/6 ~ 5.45. 

5. A value of .^ : .^ = ^ 0.28. 

6. A matrix P(A) : 


-2iX 

1 + e*'*' + e“*^ 

g2iA _j_ g, 

-2iX 

1 + e*'*' + e“*^ 

g2iA 


1 + e*'*' + 6“*'*' 

g2iA 


P(A) = 


7. An equation for greatest eigenvalue ; 2 (A) of matrix P(A) : 

— A^ 2;^(1 + 2 cos A + 2cos2A) — A 2 :(l + 2 cos A) + 1 + 2 cos A = 0. 

8. A value of deterministic diffusion coefficient: 

31A - 16 

D = —--- ^ 0.89. 

36A - 24 

9. Values of densities aj on subintervals 4,j : ai = as = ~ 2.22, 


Example 8. 1. Markov subpartition of m = 3 subintervals: 

Ik,i = [k- ^,k- 0,4,2 = [k - ^,k + ^),Ik,3 = [k + ^,k + l),k e Z. 

2. A system of two equations for ^ and the value of slope A : 


2.22, 0:2 — 


A« = |. Ai = 2 + e 

3. An equation for A : A^ — 4A — 3 = 0. 

4. A value of A : A = 2 + \/7 zz 4.65. 

5. A value of ^ ^ 0.32. 

6. A matrix P(A) : 


P(A) = 


1 + 
1 + 
1 + 


7. An equation for greatest eigenvalue ; 2 (A) of matrix P(A) : 

—A^z^ + A^ 2:^(1 + 2 cos A) + A; 2(1 + 2cos2A(l + 2 cos A)) + 1 + 2 cos A = 0. 

8 . A value of deterministic diffusion coefficient: 

^ 9A + 2 479 + 107^7 ^ 

13A - 9 894 

9. Values of densities aj on subintervals Ikj '■ ai = 0:3 = A + 0:2 = 1 + 

Example 9. 1. Markov subpartition of m = 4 subintervals: 

7fc,i = [k- ^,k- 0, 4,2 = [k-^, k),Ik,3 = [k,k + ^), Ik,4 = [k + ^,k + l),k E Z. 
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2. A system of two equations for ^ and the value of slope A : 


Ae = l, 

3. An equation for A : A^ — 2A — 2 = 0. 

4. A value of A : A = 1 + \/3 ~ 2.73. 


5. A value of = 

6. A matrix P(A) : 


_ v^-l 


0.37. 


A- = 1 + £. 
2 ^ 


/ 0 


P(A) = 


A 


^—iX 


V 


0 

0 


1 

1 

^—iX 

^—iX 


^iX 

^iX 

1 

1 


0 \ 

0 

^iX 


0 / 


7. A value of greatest eigenvalue 2 :(A) of matrix P(A) ; 
2 :(A) = ;^(1 + Vl + 2 cos A). 

8. A value of deterministic diffusion coefficient: D = 

9. Values of densities aj on subintervals 
1.08. 


0 . 11 . 


_ 3+v^ 

CXi — C14 — —g— 


0.79, «2 = as = 


Example 10. 1. Markov subpartition of m = 4 subintervals: 

4,1 = [k- ^,k- 0,4,2 = [A; - 4 A;), 4,3 = [A:, A; + 0, 4,4 = [A; + 4 A; + A), A; G Z. 
2. A system of two equations for ^ and the valne of slope A : 


Ae = l, 

3. An eqnation for A : A^ — 4A + 2 = 0. 

4. A value of A : A = 2 + \/2 ^ 3.41. 

5. A value of 

6. A matrix P(A) : 


A- = 2-A 
2 ^ 


0.29. 



/ 

e 

1 

giA 

g2iA 

1 


g-iA 

1 

giA 

0 

A 


0 

g-iA 

1 

giA 


V 

g-iA 

g-iA 

1 



9. Values of densities aj on snbintervals 4,j 
0.29. 




P(A) = 


7. A valne of greatest eigenvalne 2 :(A) of matrix P(A) : 
2 :(A) = ;^(1 + cos A + a/ 2 — sin^ A). 

8. A value of deterministic diffusion coefficient: P = / 


ai = a4 = ^ ^ 


0.21, a2 — as — 


- - 2 -A 2 

2 


Example 11. 1. Markov subpartition of m = 3 snbintervals: 

4,1 = [A^ ~ / A; — ,^ 2)5 4,2 = [A^ ~ 4) A: — ■Ci)) 4,3 = [A^ ~ 4) A: + .^i), 4,4 = [A^ + 4) A: + 4); 4,5 
[k ^21 k ^),k G Z. 

2. A system of eqnations for ^ and the valne of slope A : 


A4 = 


2 ’ 


A6 = 2-6, 


A5 = 2 + &. 
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A3 _ 4A2 _ 4A + 3 = 0. 


3. An equation for A : 

4. A value of A : A ^ 4.75. 

5. A value of ^ ^2 = 

6. A matrix P{X) : 


P(A) 


/ 

0 

0 

1 

+ 

giA 

+ 

e" 

-iX 

g2jA 

0 \ 


g-2jA 

0 

1 

+ 

giA 

+ 

e" 

-iX 

g2iA 

0 



0 

1 

+ 

giA 

+ 

e" 

-iX 

0 

g2*A 


0 

g-2*A 

1 

+ 

giA 

+ 

e" 

-iX 

0 

g2iA 

V 

0 

g-2*A 

1 

+ 

giA 

+ 

e" 

-iX 

0 

0 / 


7. An equation for greatest eigenvalue 2 :(A) of matrix P(A) : 

+ 2 cos A) — A^ 2:^(1 + 2cosA(l + 2 cos A)) — A^ 2 ;^(l + 2 cos A + 2cos2A) — z + 


1 + 2 cos A = 0. 

8. A value of deterministic diffusion coefficient: 


9. 

0:3 


Values of densities aj on subintervals Ikj : 

_ 5A^-llA+6 
16(9-A) • 


D 

tti 


81A2 + 69A - 55 

=-0.64. 

131A2 + 99A - 93 

= «5 = 4A(9-A) > a2 = a4 = 


A^+4A-3 
16(9-A) ’ 


Example 12. 1. Markov subpartition of m = 6 subintervals: 

4,1 = [k - \,k - 6), 4,2 = [k - i2,k - 6), 4,3 = [k - 6,^), 4,4 = [k,k + 6), 4,5 = 
[k + ^i,k + ^ 2 ) ) 4,6 [k + ^2, k -h Xj, k G Z. 

2. A system of equations for ^ and the value of slope A : 


^4 — 4) A,^2 —2 —4) A- —2 + ,^!. 

3. An equation for A : A^ — 4A^ + A — 8 = 0. 

4. A value A : A 4.22. 

5. A value ^ : 4 = 4 = a%- 

6. A matrix P(A) : 


P(A) 



( ° 

1 + e-*^ 

0 

0 

giA g2iA 

0 \ 


0 

g-iA 

1 

0 

giA _|_ g2iA 

0 

1 

g-2jA 

g-iA 

1 

0 

giA 

g2iA 

A 

g-2iA 

g-^A 

0 

1 

giA 

g2iA 


0 

e-i\ ^ 

0 

1 

giA 

0 


V 0 

g-iA _|_ g-2iA 

0 

0 

1 + e*^ 

0 / 


7. An equation for greatest eigenvalue 2 :(A) of matrix P(A) : 

—A'^z'^ + A^ 2;^(1 + 2 cos A) + A^ 2 :^(l + 2 cos A) + z{l + 2 cos A + 2cos2A + 2cos3A) + 2 + 


2 cos 2A + 4 cos A cos 2A = 0. 

8. A value of deterministic diffusion coefficient: 


9. Values of densities aj on subintervals I, 

^ ^ A^+A+5 

tt 3 CI4 3A2_8A+1- 




D 

Ci\ 


5A2 + 26A + 22 

=-~ 0 49 

18A2 + ISA + 56 

_ _ A^ + l _ _ 

- tte - 3A2_8 a+i’ «2 - as - 


A^+2 

3A2-8A+1’ 


Example 13. 1. Markov subpartition of m = 7 subintervals: 

4,1 = [k-\,k- 4), 4,2 = [k-iz,k- 4), 4,3 = [/c - 4, ^ - 4), 4,4 = [A; - 4 , /c + 4), 4,5 = 
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[/c + ^1, fc + ^2),Ik,6 — [^ + '^25^ + ■Cs); Ik,7 — [A: + ^ 3 , /c + |) /c G Z. 

2. A system of equations for ^ and the value of slope A : 


A6 = 6, A6 = e3, ^^3 = ^, Ai = 2-ei. 

3. An equation for A : A^ — 4A^ + 1 = 0. 

4. A value of A : A 3.968. 

5. A value of ^ = ^, ^ ^3 = 

6 . A matrix P(A) : 



/ 


g-iA 


1 

0 

0 

0 

0 

giA _J_ g2iA 

\ 




g-iA 


0 

1 

0 

0 

0 

giA _j_ g2iA 


1 



g-iA 


0 

0 

1 

0 

0 

giA _j_ g2iA 


PW = ^ 



g-iA 


0 

0 

1 

0 

0 

giA 



e" 

+ e- 

-2i\ 

0 

0 

1 

0 

0 

giA 




e" 

-t\ g- 

-2i\ 

0 

0 

0 

1 

0 

giA 



V 

e" 

-t\ g- 

-2i\ 

0 

0 

0 

0 

1 

g*A 

/ 

7. An equation for greatest eigenvalue z 

(A) 

of 

matrix 

P(A) 



AV-AV(2 + 2cosA) + l = 


0 . 










8 . A value of deterministic diffusion coefficient 

: D 

= 

4(A 

-3) 

^ 0.26. 



9. Values of densities aj on subintervals Ik,j '■ ai = aj = ^2 = «6 = f > «3 = «5 = 

(A2_3A_3)A _ _ 3 _ 5 

4(A-3) ) “4 4 4 2(A-3) ' 

The results of theorem [3] and above-considered examples allow us to give the following 
dehnition of deterministic diffusion of LDS 0 - 0 . 


Definition 6 . We say that the LDS ([I])-® has a deterministic diffusion if, for any initial 

probability measure po with bounded density, there exist a sequence of numbers cr^ > 0 

1/2 

and and 1-periodic function a{x) > 0, f a(x) dx = 1, such that the sequence of mea- 

- 1/2 

sures Hn = obtained from the initial measure by the n-fold action of the LDS, is 

asymptotically equivalent, as n —)■ cxo to a sequence of measures with densities 


Pn{x) 





2(7 


T 

71 


2 


Remark 2. The LDS ([I])-® on the entire axis generates the associated DS Xn+i = f{xn) on 
the finite segment [—A, 1 ), where the value of the function f{xn) is equal to the fractional 
part of f{x), i.e., f{xn) = {f{x)), and the function / maps the segment [—|, |) into itself. 
We denote this DS, which is a compactihcation of the LDS, by CLDS. The function a{x) > 0 
in Dehnition in] is the density of the invariant probability measure for the CLDS. 


5 Deterministic Diffusion in the Case of Transporta¬ 
tion in a Billiard Channel 

The results of numerical experiments carried out in [2] demonstrate that the deterministic 
diffusion in long billiard channels is anomalous. There are different theoretical models of 


17 










this type of diffusion transport in long channels. One of these models can be found in [T]. 
The essence of this model can be described as follows: 

Consider a long billiard channel with boundaries in the form of periodically repeated arcs 
slightly distorting the straight lines of boundaries. Assume that the upper boundary of the 
channel is symmetric to the lower boundary. Then the trajectory of motion of a billiard ball 
regarded as a material point obeying the ideal law of reflections from both boundaries can be 
studied in a channel of half width with symmetric reflections of the trajectory in the upper 
part of the channel into its lower part relative to the straight middle line of the channel. 
Hence, the reflections from this middle line can be regarded as ideal. 

We approximately assume that the reflection from the lower part of the distorted bound¬ 
ary is realized in its linear approximation. However, the normal v to the surface of actual 
reflection does not coincide with normal to the rectilinear approximation of the channel 
and is described by a known function periodic along the axis of the channel. Assume that 
Xn-i, Xni and Xn+i are the abscissas of points of three consecutive reflections of the billiard 
ball and that the vector z/(a:„) makes an angle a{xn) with the normal to the surface (see 
Fig. 2). Then the analysis of the ideal reflection at the point Xn leads to the equality of the 
angle of incidence Lpn and the angle of reflection relative to the vector u. It is clear that 
~ tan((p — a) and = tan('0 -|- a). Hence, 

Xn+i - Xn ^ ^ tan(2a) 

h 1 - tan(2a) • ' 

For large h, we get the following approximate relation from fl2^ : 

Xn+l - Xn = Xn- Xn-1 + h tail{2a{Xn)) ■ 



Figure 2: Reflections of the billiard ball in a channel of half width. 

Let the function f{x) = tan(2ax) be periodic in x with period 1. Then we arrive at the 
following model of a trajectory in the billiard channel: 

Xn+l-Xn=Xn-Xn-l+ f{Xn). (24) 

The dynamical system fl2T)l is a two-dimensional dynamical system discrete in time which 
generalizes the dynamical system ([T]). The initial conditions for xq and xi are assumed to 
be given. Let Xq = 0 and let Xi E Ik = [k — k + ^). Equation fl24)) can be represented in 
the equivalent form as follows: 

n 

Xn+i=xi-n + '^{n + l-k)f{xk). (25) 

k=l 
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Assume that the function f{x) = A-{x) in Eq. fl2^ is linear with respect to the fractional 
part of the argument. Then, for A > 1 and Xi uniformly distributed over the interval the 
terms in can be regarded uniformly distributed independent random variables. This 
leads to the normal distribution of Xn+i with variance equal to the sum of variances of all 
terms in fl25|l . 


Hence, the variance of the distribution of Xn+i is equal to 



+ 


A^ n{n + l)(2n + 1) 
12 6 


and nonlinearly depends on n oo. The deterministic diffusion in 


this billiard strip is anomalous. Thus, there are two factors leading to the anomaly of the 
deterministic diffusion, namely, the quadratic dependence of the variance on distribution of 
the initial values even for the ideal billiard and the growth of coefficients of the terms in sum 
leading to the cubic dependence of the variance of Xn+i on n. 

We now make an important remark. We have studied the distribution of positions of the 
billiard ball after n ^ oo reflections. From the physical point of view, it is more important 
to get the distribution of the abscissas of billiard balls for large values of time t because, 
for the the same period of time, the billiard ball makes different numbers of reflections for 
different trajectories. 

In the ordinary billiards, the time between two successive collisions is proportional to the 
covered distance. However, in the model of bouncing ball [13] over an irregular surface, the 
time between two consecutive collisions is maximal if the points of consecutive reflections 
coincide. The problem of investigation of the anomalous deterministic diffusion in billiard 
channels as f —)■ cxd is of signihcant independent interest. 

Note that the Gaussian density in the case of ordinary diffusion is the Green function of 
the Gauchy problem for the partial differential equation 


du 

'm 


— DAu 


0 . 


In the investigation of anomalous diffusion, we consider a differential equation with frac¬ 
tional derivatives. Equations of this kind are studied in numerous works (see m and the 
references therein). 
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